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Abstract 

'd , 

\ The charge and entropy currents across a quantum point contact is expanded as a series in 

powers of the applied bias voltage and the temperature difference. After that, the expansions of 



the Seebeck voltage in temperature difference and the Peltier heat in current are obtained. With a 
suitable choice of the average temperature and chemical potential, the lowest order nonlinear term 



in both cases appear to be of third order. The behavior of the third-order coefficients in both cases 
O ' are then investigated for different contact parameters. 
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I. INTRODUCTION 



Various aspects of the ballistic electron transport across quantum point contacts are 
studied extensively in the past. The most striking feature of this transport is the quantization 
of conductance^ at integer multiples of the conductance quantum 2e 2 /h. This phenomenon 
is usually treated with the Landauer-Biittiker formalism 3 -* 4 - which provides a transparent 
explanation for the effect. Electrons in each sub-band corresponding to the transverse modes 
in the contact contribute one quantum to the conductance if the sub-band is sufficiently 
populated. As the size of the constriction is changed by varying the negative voltage on 
split gates, which are used to define the contact on a two-dimensional electron gas, the 
conductance changes in smooth steps from one conductance quantum into the other. It is 
observed that the linear Seebeck and Peltier coefficients for these structures display quantum 
oscillation a 5 i 6 i 7 i 8 i 9 with peaks coincident with the conductance steps. 

Nonlinear transport in these systems has also been studied extensively both 
theoreticall y 10111112113114115116 and experimentall y 1 l18 Since Onsager's reciprocity relations 
connecting the Seebeck and Peltier transport coefficients loose its meaning in this regime, 
these two effects show distinctively different behavior. New peaks appear in the differential 
Peltier coefficient as the driving voltage is increased , 15116 while the thermopower does not 
change much even for very large temperature differences.— 

A major theoretical difficulty in the nonlinear regime is, due to the small size of these 
systems, finite voltage differences create large changes in the distribution of electrons around 
the contact. As a result, more involved calculations are necessary for describing the electron 
transport . 19 i 20 However, it is of some interest to analyze the nonlinear transport properties 
without taking such changes into account. The purpose of this article is to investigate the 
nonlinearities in not so commonly studied Seebeck and Peltier effects, assuming that the 
contact potential is not changed apart from the uniform shift caused by the gate voltage. It 
is hoped that this will clarify the importance of the effects mentioned above. In the following 
section, the charge and heat currents are expanded as a series in powers of the potential and 
temperature differences. Appropriate expansions for the Seebeck and Peltier phenomena are 
obtained and the series coefficients are investigated in sections 2 and 3, respectively. Finally, 
the results are summarized and discussed. 



2 



II. THEORY 



In the following we consider two electron gases connected by a quantum point contact. 
The chemical potentials and \x R and the temperatures 9l and 6 R of the left (L) and right 
(R) reservoirs are the parameters that define the whole system. The difference between the 
chemical potentials, A/x = u^ — ur, is equal to (— e)V where V is interpreted as the electrical 
potential difference between L and R. A difference in temperatures A8 = 9l — 9r as well as 
a potential difference cause electron transport which can carry both charge and heat across 
the contact. The average currents on the contact are completely determined by the sum 

T{E) = YT n {E) , 



n 



where T n (E) is the transmission probability of an electron with energy E incident from the 
nth mode. The charge and entropy currents from L to R can then be expressed as 21 

I = 2^ J°° dE(f(x L ) - f(x R ))T(E) , (1) 

Is = 2-^ J dE(s(x L ) - s(x R ))T(E) , (2) 

where 

fix) = , (3) 

s(x) = -/(x)log/(x)-(l-/(x))log(l-/(x)), (4) 

xl,r = -r— ft , (5) 

and the spin degeneracy factor is added for both currents. 

For the case of weak nonlinearities, it is useful to expand the currents in terms of the 
driving temperature and potential differences A# and V. In order to do this the variable 
of integration is changed from energy £ to a dimensionless variable denoted by x, which is 
defined as the arithmetic average of xl and x R . 

x = -{xl + xr) ■ 
This leads us to define average temperature and chemical potentials by 

k B 9 

9L~+eR~ ' (7) 

_ Qr^l + 9 L fi R 

11 = o K + o L ■ (8) 
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Here, 9 is the harmonic average of the temperatures of the two electron gases and ~p is an 
average of chemical potentials weighted by inverse temperatures. These two quantities will 
be considered as the fundamental parameters describing the contact. In other words all of 
the transport coefficients are considered as functions of these average quantities. 

With these definitions the energy variable can expressed as E = ~p + xk B 9 and the 
difference of the dimensionless x parameter is 

Afi + xk B A9 

Ax = x L - x R = — (9) 

K.BVA 

where 6 a is the arithmetic average of the temperatures on both sides of the contact 

9a = \{9l + 9 r ) . 

Finally, dimensionless driving forces are defined as 

A9 . . 

c = J~ a , (10) 

The obvious advantage of these definitions is the elimination of some terms in the power 
series expansion of the integrands in equations (JIJ and (J2J). We have 

,(-e)^ k B 9 



I = 2 ^ V 

h ^2 2m (2m+l)! 

m=0 v 1 



X 



J dx (5 + xe) 2m+1 f {2m+1) (x)T(jl+xk B 9) , (12) 
t r> k B sr^ k B 9 

Iq = 2 > X 

h ^ 2 2m (2m + 1)! 
x J dx (5 + xe) 2m+1 s {2m+1 \x)T(jL + xk B 9) , (13) 

where even order derivatives of the functions f(x) and s(x) have disappeared. This is the 
primary reason for defining the averages in Eqs. (|ZJ) and (jHJ) in this particular way. Defining 
the parameters 

f m , P = f m ,M~8) = dxx p f m \x)T(jl + xk B 9) , (14) 

which are only functions of the contact parameters /Z and 9, the currents can be expressed 
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as 

oo 2m+l 



m=0 p=0 ^ v ^' 



is = 2 ^k B eyy 1 [ -^^ . (i6) 

m=0 p=0 ^ v ^ 7 

This is the desired expansion of currents in terms of the driving forces e and 5 with the 
coefficients being functions of the average quantities ft and 9. 

One notable property of the equations (fT5j) and (fTfij) is that only the odd powers of the 
driving forces combined together appear in those expressions. This implies that if both 
driving forces change sign e —>■ — e and 5 — > —5 then the charge and entropy currents change 
direction. Including only up to the third order terms in the expansions we have 

/ = 2^-k B 9(f 10 6 + f n e 
h 

+ Y A (/3o5 3 + 3/ 31 5 2 e + 3/ 32 5e 2 + / 33 e 3 ) + ■ ■ ■ ) (17) 

I S = 2^- (f n 5 + foe + ^ ((/ 31 - 2f 20 )5 3 + 3(/ 32 - 2f 21 )5 2 e 

+3(/ 33 - 2f 22 )6e 2 + (/ 34 - 2 / 23 )e 3 ) + • • • ) (18) 

These equations give the currents for arbitrary values of the temperature and potential 
differences. However, measurements are rarely carried out for arbitrary A9 and V. Electrical 
conductance and Peltier effect measurements are carried out at isothermal conditions while 
the thermal conductance and Seebeck effect measurements are done with zero electrical 
current. But, the equations above is a starting point for each particular phenomenon. In 
the following, only the Seebeck and Peltier effects are investigated. 



III. SEEBECK EFFECT 



In the Seebeck effect, a temperature difference creates a potential difference across the 
point contact when there is no electrical current (1 = 0). This potential difference can be 
expressed in dimensionless form as 

-6 = cj ie + cT 3 e 3 + a 5 e 5 + ■■■ (19) 
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where the first two coefficients are 



hi 



(20) 




1 



(/33 - 3/32^1 + 3/31 cti - f 30 af 



) 



(21) 



In terms of V and A9 the series expansion is 



V = StAO + S 3 A9 3 + S 5 A9 5 



+ ... 



(22) 



where 



5, 



A; B 1 



m = 1, 3, 5, . . . 



m 




Appearance of only the third order terms in Eq. ()22|) implies that when the temperatures of 
the two reservoirs are exchanged (in other words the sign of A9 is changed without changing 
9 a and 9), the induced potential difference due to the Seebeck effect is reversed. 
The nonlinear terms in Eq. ()22j) becomes significant when 



It is possible to get a theoretical estimate of this quantity in the small temperature limit, 
when k B 9 <C El, where El is the energy range where T(E) changes by one. In this case, 
the Taylor series expansion 

T(jl + xk B 9) ss T(JI) + xk B 9T'{jl) 
in Eq. (jl4j) gives the following approximate expressions for <j\ and 03 

7r 2 r - 7r 2 r' - 

° x ~ !sT kBe ' as ~ 12 ~r kB6 ' 

The threshold level for nonlinearity is then 



Since A# can never go above this level, the nonlinearities in the Seebeck effect are always 
smalL 1§ For this reason, the expansion (}2*2*|) is appropriate for almost all nonlinear cases. 
For the opposite, high temperature limit, numerical calculations of the Seebeck coefficients 
indicates that the threshold expression given above does not change much. 




A$threshold ~ 26 A = &l + 6r 
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FIG. 1: (a) The third order Seebeck coefficient, <r 3 = {—e)9\/kB S3, is plotted as a function of 
average chemical potential JI for Lv y /io x = 6 and ksO/fv^y = 0.01, 0.02, 0.04, 0.08, 0.105 and 0.125 
(from bottom to top). Each curve is shifted by 0.05 units for clarity. (b)For comparison the linear 
Seebeck coefficient, o\ = (—e)/kB Si, is plotted for the same set of parameters. Each curve is 
shifted by 0.02 units and the temperature increases from bottom to top. 



As for the general behavior of 5*3, we calculate it for a contact denned by the saddle 
potential 

\rl \ 1 2 2 1 1 2 2 

V{x,y) = ~-mu x x + -muj y y . 
For this case the energy dependent transmission probability for the nth transverse mode 



[n = 0,1,2,...) is 



T n {E) 



1 



l + «p(-&^-^y("+5)]) 
In Fig. [H 5*3 is plotted against Ji for this potential. At sufficiently low temperatures, third 
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FIG. 2: (a) The third order Seebeck coefficient, 03 = (—e)0\/kB S3, is plotted as a function 
of average chemical potential fx for ksO/hujy = 0.04 and uj y /uj x = 1.5, 3, 6 and 12 (from top to 
bottom) respectively. Each plot is shifted by 0.025 units for clarity. (b)For comparison, the linear 
Seebeck coefficient, a\ = {—e)/kB S\, is plotted for the same set of parameters. Each plot is shifted 
by 0.2 units and uj y /uj x ratio increases from top to bottom. 

order Seebeck coefficient, S3, has single peaks coincident with the peaks of Si. When the 
temperature is increased, these peaks start to split into two. This change happens around 
ksO/fkJx ~ 0.08. It is observed that the distance between the peaks is proportional to the 
temperature. For this reason, with increasing temperature, the structure develops into two 
separate peaks. Also, the widths of the peaks increase proportionally with the temperature. 
Inevitably, when the temperature is increased further (around ksO/huy ~ 0.08), each peak 
of the pair starts overlapping with the peaks of the neighboring steps. For this reason, in 
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this high temperature regime the nonlinearity in the Seebeck effect becomes more significant 
away from the steps (at the plateaus of the electrical conductance). Same graphs are shown 
in Fig. El for different values of w y /uj x ratio. It can be seen that S3 has single peaks for small 
values of LO y /u x ratio (around lo v /u x ~ 1), and peak splitting occurs for larger values of the 
ijjyjijj x ratio. 

In all cases it can be seen that S3 is always negative (03 is always positive) and never 
changes sign. It implies that the nonlinearity increases the generated Seebeck voltage further 
than the linear term alone suggests. Note that this feature of S3 is not apparent from its 
definition, Eqn. (|2Tj) . This appears to be a model dependent feature. Especially if T(E) 
may decrease for some energies, S3 may display sign changes. But for the saddle potential 
model and for all parameter ranges investigated in this study, S3 is found to have the same 
sign. 



IV. PELTIER EFFECT 

The Peltier heat is defined as the heat carried Q = 61$ by the charge current I at 
isothermal conditions (9l = Or = 0). The expansion of the Peltier heat and the charge 
current in terms of the 5 parameter is 



Q = 2^- [f n 5 + -(/31 - 2fa)5* + — (/ 5 i - 4/ 40 )5 5 + • • • j , (23) 
/ = 2^(/ 10 * + l/^ + i l 5 /^ + ...) - (24) 
Both of these expressions can be used to expand Q as a power series in the current / 

q = n a i + n 3 / 3 + n 5 / 5 + • • • , (25) 

where the first two terms of the expansion are 

jj k B 9 fu 

= ' (26) 

■pr h /lp(/31 ~ 2/20) ~ /ll/30 /r, 7 N 

113 " {-efk B 9 96A 4 ' KM) 

The appearance of only the odd powers of the current in the expansion of Q signifies the 
reversible character of the Peltier heat. The coefficient Hi is for the linear Peltier effect, 
which is related to Si through the Thomson- Onsager relation by II i = 9Si. 



The plots of n 3 are shown in Fig. |3] and |U for the saddle potential model as a function of 
JL for different values of parameters ksO/huy and u y /u x , respectively. For low temperatures, 
II3 is non-zero only around the steps of the conductance. But, in contrast to 5*3, it displays 
a change of sign for all parameter values. In particular II3 has opposite sign at the peaks 
of 111 = OS. This behavior is an indication of the peak splittin g 1 ^ 1 ^ behavior of the Peltier 
coefficient under nonlinear currents. In other words, with nonlinear currents, the Peltier 
heat decreases at the peaks of the linear Peltier coefficient, but increases at the foothills of 
these peaks. Similar to S3, II3 is extremely small at the plateaus of the conductance for 
small temperatures, but when the temperature is higher (comparable to frw y ) it also becomes 
significant at the plateau region. Finally, II3 is significant only around the first few steps. 
At higher steps, it is observed that the peak heights are inversely proportional to the cube 



To estimate the threshold level for nonlinearity, we use the following approximations valid 
in small temperature limit 



in Eq. (j2HJ). Therefore, the nonlinearity sets in when the driving potential difference is of 
the order of eVthreshoid ~ El- Since it is possible that the driving potential difference on 
the contact can easily exceed this threshold level, in these highly nonlinear cases it will not 
be reasonable to use only a few terms of the expansion in Eq. (|25|). However, for weakly 
nonlinear cases, the expansion above might be useful. 

High-order nonlinearity in Peltier effect at small temperatures 

As it was discussed above, highly nonlinear cases cannot be treated appropriately by the 
power series expansion discussed here. For this case, we need to have a better method for 
evaluating the heat and charge currents passing through the contact. We consider only the 
isothermal case appropriate for the Peltier effect. The charge and entropy currents for this 



QfT(/Z). 
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FIG. 3: The third-order Peltier coefficient II3 (in arbitrary units) is plotted as a function of 
average chemical potential /I for uj y /uj x = 6 and different values of temperatures {k B /huj y values 
are indicated in the figure). 
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FIG. 4: The third order Peltier coefficient H3 is plotted as a function of average chemical potential 
/Z for k B 8 jhjjy = 0.04 and different values of u> y /uj x whose values are indicated in the figure. 



case can be expressed as 

I = T— L \ dx(-f'(x))[A(pL L + xk B 9)-A(fjL R + xk B 9)} , (28) 
J -00 
k f°° 

Is = 2-^ / dx(-xf'(x))[A(pi L + xk B 6)-A(pi R + xk B 9)} , (29) 
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where A(E) is the energy integral of T(E), 



A(E) = / T(E)dE 



Assuming small temperatures (ksO <C El), the integrands can be expanded as 



AQj, + xk B 9) Rj A(ji) + xk B 9T(fi) 



Keeping only the lowest order terms the currents can be expressed as 




(31) 



(30) 



As was discussed by Bogachek et a/. , 15 i 16 the differential Peltier coefficient can be expressed 
as (assuming constant ~p) 



The peak splitting effect of the nonlinearity can be seen from this expression. When the 
potential difference across the contact is less than El, the individual peaks of T'(/j,l) and 
T'(fiR) will join in a single peak observed in the linear Peltier effect. However, if the po- 
tential difference is more than El, the contribution of these two terms can be distinguished 
since they will form two separate peaks. The distance between the peaks, then, will be 
proportional to the applied potential difference. 

V. CONCLUSIONS 

The expansions of the charge and entropy currents as a power series in temperature 
and potential differences are obtained, assuming that the transmission probabilities are 
unchanged by the nonlinearities. The main advantage of this particular expansion is, through 
a different definition of average chemical potential, ~p, and temperature, 9, some particular 
terms disappear from the expressions. The Seebeck and Peltier effects are investigated as 
special cases and it is found that the lowest order nonlinearities are of third order in both 
cases. 

In the case of the Seebeck effect, S3 is found to have the same sign as S\. Although 
at low temperatures S3 is found to be simply proportional to Si, its peaks split into two 
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at high temperatures. If ksO is comparable to the energy difference between the successive 
sub-bands, these peaks may join with the peaks of the neighboring steps, creating an unusual 
appearance where S3 has maxima at the plateaus of the conductance and minima at the 
steps. In all cases, it is found that the nonlinear signal is small compared to the linear one. 

For the case of the Peltier effect, II3 changes sign as the gate voltage is changed for all 
parameter values. The main shortcoming of the expansion developed here is that in this 
case the potential difference driving the current may be chosen above the threshold level 
for nonlinearity. In such a case, the expansion is useless as more and more terms have to 
be added up to obtain the correct response. In the small temperature limit, an alternative 
expression has been developed for the differential Peltier coefficient that is also valid for 
highly nonlinear cases. 
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